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^ ' We consider an array of XX spin- 1/2 chains coupled to acoustic phonons and placed in a magnetic 

field. Treating the phonons in the mean field approximation, we show that this system presents a 
first order transition as a function of the magnetic field between a partially magnetized distorted 
state and the fully polarized undistorted state at low temperature. This behavior results from the 
magnetostriction of the coupled chain system. A dip in the elastic constant of the material near the 
. saturation field along with an anomaly in the magnetic susceptibility is predicted. We also predict 

O ' the contraction of the material as the magnetic field is reduced (positive magnetostriction) and the 

ecu ■ 
reciprocal effect i.e. a decrease of magnetization under applied pressure. At higher temperature, 
the first order transition is replaced by a crossover. However, the anomalies in the susceptibilities 
in the system near the saturation field are still present. We discuss the relevance of our analysis in 
"j^ ' relation to recent experiments on spin-1/2 chain and ladder materials in high magnetic fields. 
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The subject of Bose-Einstein Condensation (BEC) has become tremendously active both experimentally and theo- 
' ^ ' retically in the recent years due to the observation of BEC in trapped atomic gasesi Another way of obtaining a BEC, 
predicted theoretically some time ago, is to place a quasi one dimensional spin gap systems under a strong magnetic 
£gj(^j2j3j4j5j6j7j8j9 causing the formation of a Luttinger liquid of magnons. This Luttinger liquid possesses a quasi-long 
range order, which is a one dimensional precursor of the BEC^Siiiii^ In the presence of any three dimensional coupling, 
this quasi-condensed state develops a long-range order, which corresponds to a BEC of the magnonsii^iM This conden- 
sation has been observed experimentally^'^ in spin-1 chain materials, and more recentljiiSiiLi^ in Cu2(C5Hi2N2)2Cl4 
I (CuHpCl), a material originally believed on the basis of its thermodynamics to be formed of weakly coupled two 
J — leg ladders. However, it is now known that CuHpCl is in fact formed of dimers coupled in a three-dimensional 
I network^^, and that the Bose condensation is a condensation of a two or three-dimensional magnon Bose fluidi^S 
' Recently, such a Bose-Einstein condensation of magnons has also been observed in other coupled dimers systems 
^""^ \ such as TlCuCls^^, KCuClgi^^, Cs3Cr2Brg^ and BaCu2Si206-^'' Theoretical description taking into account the two- 
i dimensional nature of these gapped systems have been develop edSSiSLSLSS in recent years. They are based on the 
bond operator theory (BOT) description of the coupled dimer systemi22i2Si2i The interactions between the condensed 
triplet bosons are treated by the Hartree-Fock-Popov (HFP) approximation^^, and a rigid lattice is assumed. However, 
more recently, some experimental evidence for a lattice distortion associated with the Bose Einstein Condensation 
has been obtained>Si2^i2^ Older reports of specific heat anomalies also pointed out to a lattice distortion associated 
1^ ■ with the transitioni2& It has thus been suggested that the distortion observed in experiments could be the generalized 
spin-Peierls transition predicted originally for spin ladder systems under field^ However, similarly to the conven- 
O tional spin-Peierls transition^SiSS, in a generalized spin-Peierls transition the lattice distortion results from an optical 
J-^ phonon becoming static. As a result, in the generalized spin-Peierls transition, the lattice constants of the crystal do 
^ t not change, but superlattice peaks are visible in X-Ray diffraction, in contrast with experiments. '^^'^^ In fact, in the 
case of a condensate of magnetic excitations coupled to a deformable lattice, another instability than the generalized 
spin-Peierls transition is possible, namely magnetostriction. Although magnetostriction has been mostly discussed 
in the context of ferromagnetism^Sf^ it can also be observed in paramagnetic or antiferromagnetic systemsi^ In the 
case of magnetostriction, the lattice parameters change as the system becomes magnetized, but the number of atoms 
per unit cell does not change, which means that no new reflection appear in X-Ray diffraction experimentsi^^ In 
a ferromagnet, magnetostriction is known to turn the second order Curie transition into a first order transitions^ 
Thus, the magnetostrictive effects may also explain the first order transition as a function of magnetic field observed 
in magnon Bose-Einstein condensation experimentsi2i^ii2^ Besides the Bose-Einstein condensation of magnons in 
spin gap systems, magnetostriction effects are also relevant near the saturation field in the spin-1/2 chain material 
Cu(II)-2,5bis(pyrazol-l-yl)-l,4-dihydroxybenzene (CuCCP)44 In fact, the two effects can be closely related as a ladder 
system under strong magnetic field can be mapped for strong rung coupling on an anisotropic spin-1/2 chain. In 
the present paper, we will discuss a model of magnetostriction in ID XX spin chains coupled by 3D acoustic phonons. 
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We will study the effect of phonon coupling at the mean field level and show that the Bose Einstein Condensation 
transition of the magnons becomes a first order transition at low temperature. The origin of the first order character 
of the transition can be traced to the hardcore repulsion of the magnons, which prevents the collapse of the conden- 
sate. Treatments based on Hartree-Fock-Popov mean field theory neglect the hardcore constraint and thus fail in the 
regime where the first order transition is obtained. This stresses that a correct treatment of the hardcore constraint 
is necessary to describe the first order transition observed in experiments^ 

The plan of the paper is the following. In Sec^l we introduce a model of an array of two leg ladders coupled 
with classical phonons as a minimal model for describe the BEC in a system of hard core bosons coupled to a lattice 
deformation. From this model we deduce an effective XXZ spin chain Hamiltonian that can be rewritten in terms 
of a free fermion Hamiltonian by the Jordan- Wigner transformation. In Sec lIIII we focus on the zero temperature 
case. By a mean field theory, we discuss the first order phase transition at varying the external magnetic field and 
the pressure. In particular, we derive a Ginzburg-Landau expansion of the free energy valid for a very stiff lattice 
which allows an analytical discussion of the behavior of the magnetization and of the lattice parameter as a function 
of the magnetic field. We also calculate and discuss the behavior of various susceptibilities in the distorted case, such 
as the elastic modulus, the magnetic susceptibility and the magnetostriction coefficient that all present anomalies at 
the transition. In Sec II VI we turn to the positive temperature case and deduce first a Ginzburg-Landau expansion 
which allows us to predict the temperature above which the first order transition is replaced by a crossover. All 
the susceptibilities are then recalculated at finite temperature and some results compared to recent experiments. — 
We also consider the behavior of the dilatation coefficient and of the specific heat. Finally in Sec. we give the 
conclusions and perspectives for future work. 



II. MODEL 



In order to grasp the origin of the first order character of the BEC in a spin gap system under field coupled to 
acoustic phonons, we want to build a model that can be solved with a minimum of approximation. We will start with 
a an array of two leg ladders coupled with classical phonons. Its Hamiltonian reads: 

^ 

n,m n,m 

where Jy, J|, J± > are the exchange constants, K is the elastic constant, and is a spin-phonon coupling constant. 

We have neglected the kinetic energy of the phonons, as is done usually in treatments of magnetostriction4& In the 
case of a strong rung interaction J±, it is possible to reduce the size of the Hilbert space by considering only the 
singlet rung state and the triplet state of spin 5^ = +1 and derive an effective XXZ spin chain theoryiiiiS*^ The 
resulting effective Hamiltonian reads: 

n.m 

where the pseudospin operator is acting on the reduced Hilbert space, 5n,m = ??(itn+i,2m — Un,2m) (we are assuming 
Un+i,2m-Un,2m. = Wn+i,2m-i-i -u„,2m-i-i as wc are Only Considering magnetostriction cflFects) , K = K/if, A = J|/(2J||), 
and h' = H — J± — ^^f/^ effective magnetic fieldi^ The magnetization per rung of the ladder system is 

TO = {{S^2m + 'S'n.2m+i)) ^ {"^n.m) + 1/2- In the casc of SU{2) invariant spin ladders (J|| — J|p, one recovers 
A = 1/2 The case of CuCCP is also described by the model defined by Eq. ((2Jl with A = 1 and the magnetization 
given by (T^)- It is interesting to note that the operators = {T^ ^zLiT^ ^) satisfy the same commutation relation 
as hardcore boson creation and annihilation operators4£i^ As a result. The Hamiltonian |5J also describes a system 
of hard-core bosons with nearest neighbor repulsion and a boson-phonon interaction. The magnetic field h' plays the 
role of the chemical potential of the hardcore bosons. When h' is sufficiently negative, the density of bosons is zero. 
In the original problem described by Eq. (Q) this corresponds to the spin gap state obtained for H <^ J^. When h' is 
less negative, the density of hardcore bosons in the ground state becomes nonzero. As our system is one-dimensional, 
this leads to the formation of a quasi-Bose condensate of hard-core bosons, in which the order parameter of the Bose 
condensation possesses quasi-long range order in the ground state. Hence, the transition between the spin gap 
state and the magnetized state in coupled dimer systems^! can be viewed as the Bose condensation of hard core triplet 
bosons. In the case we are considering, we have to consider the effect of the phonons on this Bose condensation. Let 
us first show that the HFP approximation can be inadequate to deal with the magnetostriction effects resulting from 
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the spin-phonon interaction^^. Indeed, if we consider bosons with an on-site repulsion described by the Hamiltonian: 



(3) 



we notice that in one dimension and in the hniit U ^ oo, the Hamiltonian Q becomes equivalent to the Hamilto- 
nian ^ with J — t/2, V = J5, and n = h' . Let us make make the HFP approximation bi (bi) = A, Sij S in ||2Jl. 
The resulting ground state energy reads: 



E ^ N 



= N 



-zt{l + S)\X\^ - + U)\X\^ + {U + zV)\X\^ + ^5^ 



-{p + U + zt)\X\^+ [u + zV- 




1.1 + U + zt 
2(u + zV-^] 



Ct, (4) 



where z is the coordination number of the lattice, 
from below, and the minimization of Q leads to: 



For U + zV — > 0, the variational energy Eq. ^ is bounded 



^i + U + zt 
2U + 2zV- 



-eiji + U + zt), 
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zt II + U + zt 
K2U + 2zV~'^ 



-ein + u + zt), 



(5) 
(6) 



where Q{x) is the Heaviside step function. From Eqs. lO-lO, we obtain a second order phase transition as /i > 
^{HFP) _ _ jj ^ jj ^ _ z_t_ ^ ^Yyq variational energy Q becomes unbounded from below. In that 
case, the HFP approximation leads to un unphysical result |A| = oo. The origin of this unphysical behavior when 
the spin-phonon interaction becomes too strong is easily seen. The HFP approximation reduces the spin-phonon 
coupling to an effective attractive interaction between the bosons of order 1? jK, and when this attraction overcomes 
the repulsion U + zV, the ground state of the mean-field Hamiltonian becomes ill-defined. It is possible to propose 
an ad hoc modification of the HFP approximation that takes into account the hard-core constraint. In the limit 
J7 — !■ cx), one must have bjbi < 1. This constraint can be satisfied by imposing |Ap < 1 on the order parameter and 
using as the variational energy Eq. Q with U = 0. Then, one obtains that for zV < t'^/{2K), there is a first order 
transition, with A = and 5 = for ^ < -zt + zV - t^ /{2K) and A = 1 and 5 = zt/K for > -zt + zV - t^/{2K). 
The modified HFP approximation suggests that first-order magnetostriction transition caused by the condensation 
of triplet excitations can be observed also in d > 1 provided the spin-phonon interaction is strong enough. The 
preceding discussion suggests that in order to describe correctly the first order magnetostriction transition associated 
with the Bose condensation of the triplets, we need to treat exactly the hardcore interaction of the triplets. In the 
quasi-one dimensional system of Eq. 0, this can be done in principle exactly thanks to the integrability of the XXZ 
spin chain^. The mean field theory of the magnetostrictive transition then proceeds^ by writing = (5Vn, m and 
minimizing the free energy with respect to 5. The parameter 5 measures the relative change of the lattice parameter 
5 oc Aa/a. In the most general case, the free energy of the general XXZ chain in the presence of magnetostriction 
have to be obtained from the Thermodynamic Bethe Ansatz (XBA)SiiSiSi^. However, given the complexity of the 
TBA equations, it is better to consider the simplest case of A = i. e. Jif = (XX chain) as in the original theory 



of the spin-Peierls transitionSSiS. In that particular case, by the Jordan- Wigner (JW) transformatioi 



56. 



T+ = 



«n cos TT ^ 

\ 7n<n / 



T,l = aia„ - 1/2, 



(7) 
(8) 



the chain can be mapped to a chain of free fermions and the calculation of its thermodynamics is much simplified. 
We expect that this mapping is not going to affect the qualitative properties of the chain under field. In fact, when 
the magnetization per site of the XXZ chain is near saturation, i.e. when the density of triplet excitations created on 
the ladders is low, the effect of the short range interaction measured by A becomes small. This is seen in particular 
in the behavior of the Luttinger exponent K which goes to 1 in that low density limit^^. Thus, we expect that the 
magnetostriction effects that happen near the transition between the singlet state and the polarized state will be 
properly described by the model |2Jl with A = 0. However, for a quantitative approach to the full magnetization 
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curve, the use of TBA equations would become necessary. Even then, the insight gained from the study of the XX 
chain wiU be useful to devise the correct treatment of the TBA equations combined with the self-consistency condition 
on S. Returning to the XX chain, after the JW transformation, and going to Fourier space, the Hamiltonian of the 
chain reads: 

H = - J(l -5)J2 coskalak + ~ h{alak - 1/2) (9) 

k 

where we have written J = J|| and h = h' to simplify the notation. The operator ak annihilates a fermion of 
momentum k. To study the thermodynamics, we consider that the system is under fixed pressure, and work with the 
Gibbs free energji^. The variational Gibbs free energy reads: 

G(p, T, h;6)^^6'+pS+^-T r ^ ln(l + ^-((i-^)^^^)-")/^), (10) 
2 2 211 

where P is the pressure, p = rjP is the reduced pressure and we have units in which % — = ^- The value of 5 is 
obtained by minimizing G(p, T, /i; 5) with respect to 5. This yields the equation: 

/" dk e(k) 
^^ ,Kl-Sip.T.H)Mk)-WT + r ^^^^ 

and the Gibbs free energy is then Q {p, T, h) — G{p, T, h; S{p, T,h)). In the following sections, we will study the Eq. l(TT)) 
at zero temperature, and exhibit the first-order magnetostrictive transition as a function of the applied field. We will 
also discuss the effect of the magnetostrictive effects on various susceptibilities of the system. Then, we will turn to 
the effect of a positive temperature. 



III. ZERO TEMPERATURE 



In this section, we first show that at T = and within mean field theory, a first order phase transition as a 
function of the external magnetic field is obtained in the array of ID XX chain coupled to acoustic phonons. Then, 
using a Landau- Ginzburg expansion^®, valid in the case of a stiff lattice, we obtain analytically the expression of the 
discontinuity in the magnetization and in the lattice parameter. This allows us to represent the hysteresis cycle of the 
magnetization as a function of the applied magnetic field. Finally, we discuss the anomalies in the elastic constant 
and the magnetic susceptibility near the transition as well as the magnetostriction coefficient. 



A. Minimization of enthalpy at zero temperature 



1. First order transition 



In order to study the equilibrium value of the relative change in the lattice parameter S, we need an expression of 
the Gibbs free energy as a function of S valid for T = 0. At zero temperature, the Gibbs free energy l(Tn)l reduces to the 
enthalpy H{p, h; S) = Eo{h; d) +pS, where Eq is the ground state energy. The ground state of the Hamiltonian ^ is 
simply the Fermi sea of pseudofermions filled up to the chemical potential h. The Fermi wavevector of pseudofermions 
is given by —J{1 — S) cos kp = h, and the magnetization M = (T^) can be expressed as a function of kp as: 

M=^-^. (12) 
Using the magnetization M instead of the Fermi wavevector, we can rewrite the ground state of the system as: 

E{M,h,5) = --(1 - 5)cos(7rM) -|- —S^ - hM (13) 

TT 2 

The enthalpy must be minimized with respect to S for fixed h,p giving: 

KS+P+- cos(7rM) = 0, (14) 

TT 
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which is the zero temperature Hmit of Eq. Hll() , and the magnetization is given by: 



J(l-(5) simrM - h = 0. 



(15) 



Both Eqs. H14|l - p5|l are vahd only for \ AI\ < 1/2. For M — ±1/2, the kinetic energy of the pseudofermions vanishes, 
and Eq = — hM. In this saturated regime, one has obviously 6 — —p/K independent of the magnetization. 

In the regime of non-saturated magnetization, 5 varies with the magnetization and the magnetic field, implying the 
existence of magnetostrictive effects. For \M\ < 1/2, it is possible to use the equation H15|l to eliminate M and write 
the ground state energy as a function of 5 only. We then obtain the final expression of the ground state energy as: 



Eoi6, h) 



J 



h ( h 

1 ""'"'^ VW^) 



K- 



ii\h\ < J\l-S\ and: 



h/2 



(16) 



(17) 



if \h\ > J\l — 6\. The Eq. H16I) describes the XX chain not fully polarized, whereas Eq. H17|l describes the fully polarized 
XX chain, the function Eo{S,h) is plotted as a function of 6 for different values oi h/J on Fig.^ As is evident on 
Fig. n h is increased above /ic, the minimum existing for 6 = 6c < acquires a higher energy than the minimum 
at (5 = (see Figs. ^b)-(c)-(d)). As a result, at h = he, the parameter S jumps from Sc to 0. Eq. H15|) then implies 
that the magnetization also jumps from a value M = to M — 1/2. In other words, the system presents a first 
order transition at T = as a function of the magnetic field. In physical terms, the origin of this first order transition 
is easy to grasp. When the system described by is in the spin gap state, it does not contain any triplet. When 
the system starts to be polarized, the triplets first occupy the states with the most negative kinetic energy. Due to 
the hard-core repulsion, instead of being all in the lowest energy state, they form a Fermi-sea. By contracting, the 
lattice increases the exchange energy and thus makes the kinetic energy of this Fermi sea more negative. In turn, this 
favors the condensation of a larger number of triplets in the Fermi-sea, and a further decrease of the energy of the 
triplets. This process is only limited by the loss of elastic energy coming from the deformation of the lattice. The 
natural consequence is that as soon as the magnetic field overcomes the gap, a finite density of triplets appear. Due 
to the fact that the triplets form a Fermi sea, their Bose condensation is incomplete and this limitates the density 
of triplets. In the magnetized state, the system of coupled ladders possesses a quasi-long range magnetic orde r^i^i^^ 
with: 



I qZ qZ 



(-)"-" cos[7r(l -2TO)(n-n') 



|n-n'|3/2 
cos[27r77i] 



(18) 
(19) 



ifi = M-\-l/2 being the magnetization per rung of the ladder. The power-law decay of the correlations of the staggered 
transverse magnetization H18I) can be interpreted a a quasi Bose condensation^^ of the hardcore bosons of the ladder 
in the lowest energy state of momentum (7r,7r). 

In fact, even in a quasi-one dimensional array of coupled chain, a long range order of the BEC order parameter 
can develop. The reason is that in any real system, besides the spin-phonon coupling, there exists also an interladder 
exchange coupling J'Sn,2m-i ■ 82,2™+!- If this interladder coupling is small enough compared to the rung coupling 
, the singlet phase is stabls^, and no antiferromagnetic ordering is observed in the absence of an applied magnetic 
field. In the magnetized state however, this interladder coupling gives rise at T = to a long range orderingof the 
spins^'^. In our problem, the Luttinger exponent being equal to one, it is easily seen using the results of Ref. 'iS'that 
the staggered magnetization in the plane perpendicular to applied field behaves as (J'/J||)^/^. Therefore, one should 
also expect a discontinuity of the staggered magnetization from to (J'/ J)^/^ as the coupled ladder system becomes 
magnetized. This discontinuity is the signature of the BEC long range order of the triplet bosons in the magnetized 
state. 
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FIG. 1: The ground state energy as a function of S for different values of the magnetic field h. (a) For low field, there is a 
single minimum at 5 < 0. (b) For a higher field h > /ic<, a second minimum, at a higher value of the energy is present. This 
second minimum corresponds to a metastable state, (c) For a critical value of the magnetic field, the two minimum become 
degenerate, (d) Above the critical field, the absolute minimum of the ground state energy is for 5 — 0. A metastable minimum 
exists for 5 < 0. (e) For h = /ic> , the minimum at 5 < becomes unstable, (f) Above /ic> , the only minimum is for 5 = 0. In 
the range of /ic< < S < ftc> , hysteresis in S as h is varied is expected. 



2. Analytic study of the minimum of the enthalpy 



We now turn to an analytic study of the dependence of the parameter S on the magnetic field. In terms of S, the 
condition (|14|) for the existence of the minimum can be rewritten as: 



KS + p 



TT 



'J2 



This equation is rewritten as the foUowing equation of the fourth degree: 



(20) 



P 



1-4^+^ - 



K 



P 



K 



J 
ttK 



6^ + 2 



J 



P_ 
K 



K 



It can be solved by the method of Cardan as exposed in Ref. |59| p. 17, sec. 3.8.3. The dependence of S on h is 
represented on figure Inland the resuhing magnetization on figure |3| 
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It can be shown that the lattice parameter given by the solution of (|21|l has a discontinuity as a function of the 
magnetic field for a critical field hf.^ given by: 



J 



I + Zii ' ^ I' p 

8 AK 



4 V TTif / 



For a field just below ft,c> , the value of 6 is given by: 



1 + ^ 



(22) 



(23) 



and above /ic> , the value of S is simply 5 = as there are no physical solutions to the equation l|21|l . The behavior 
of the enthalpy as a function of S ior h = hc^ is represented on Fig.n}(e). 
For S = Sc^, the corresponding magnetization is M^^ such that: 



cos(7rAfc 



7r{K +p) + Tr^(if+p)2 + 8J2 



(24) 



As with all first order transition, the solution of Eq. (|21|l may only describe a metastable state and not the true 
minimum of the enthalpy, and the true thermodynamic transition happens for a field he lower than hc-^ . This is clearly 
seen on Fig.^ As a result, the field hc^ only corresponds to an extremity of the hysteresis cycle of the magnetization 
loop. The other extremity of the hysteresis cycle is obtained when the minimum at 6 — becomes unstable, i.e. when 
the second derivative of the enthalpy w.r.t. 6 becomes negative for 6 = 0. The corresponding critical field is given 
by ^c< = "^(l -^p/K). Obtaining analytically the true thermodynamic critical field is difficult in the general case. 
However, in the case of a stiff lattice this can be done easily using a Landau-Ginzburg expansion. 
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FIG. 2: Striction for T = 0, P — 0, K = 1, J = 1. This is obtained from the solution of Eq. I|21^ . It does not correspond to the 
thermodynamic striction which presents a jump for h = he < hc^ ■ 
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FIG. 3: magnetization for T = 0, P — 0, K = 1, J = 1. This is obtained from the solution of Eq. H21|l and this is not the 
thermodynamic magnetization. 



3. Landau- Ginzburg expansion at T — 



A more detailed description of the first order transition can be obtained in the limit of very small magnetization 
jump, when it is possible to expand the ground state energy in powers of M — 1/2. Expressing the ground state 
energy as a function of M only, using the Eqs. (|13|l and H14|) . we find: 



E{M) 



J2 



(ttM) - + cos(ttM) - hM 

TT V K/ 



P 

2K 



(25) 



It is convenient to substract the energy of the fully polarized system E{M = 1/2) = —h/2 — /{2K) to E{M) and 
work with the quantity: 



Eo{M) = 



J2 



(nM) + cos(7rM) 



2Tr^K 



cos 



-hi M - 



1 



Expanding Eq. (|26|l to third order in M — 1/2, we obtain: 



M--) - — IM-- 
2 2K \ 2 



M - - 
2 



(26) 



(27) 



It is important to note that since M < 1/2, the negative coefficient of the cubic term does not lead to any instability. 
Minimizing the ground state energy with respect to M, one finds: 



J 



(-1) 



- h 



M 

K \ 2 



(28) 



For h/ J > (1 -^p/K) + ,P /[2'k'^K{K +p)] ~ /ic> , this equation has no solution. In this regime, Eq{M) is a uniformly 
decreasing function of M and the minimum is obtained for AI — 1/2 corresponding to the fully saturated case. For 
h/J < {1+p/K) + J'^/[2tt^K{K +p)], the equation (gHl admits two solutions: 



M. 



</> 



1 i±^ii)-^^n^ + ^-i]i^+p/K) 

2 ^Mi + f) 



(29) 
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When h/J < 1 +p/K = hc^, the solution M> becomes larger than 1/2 and thus unphysical . In this regime, M< is 
an absolute minimum of Eq(M). Let us first consider the intermediate regime /ic< < h/J < h^^ . In this regime, there 
is a minimum of Eq{M) for M = Af<, a maximum for M = A/> and a second minimum for M = 1/2. The exchange 
of stability of these two minima is at the origin of the first order transition. We have i?o(l/2) = 0, therefore, the 
stability of the minimum at M = M< depends only on the sign of £'o (A/< ) . Using the equation H28() , we can rewrite: 



1X2 



^"("<' = r<-V L2A- ■ 3 r ■ lay- i 

The stability of the two minima at Af =1/2 and M — Af< is thus exchanged when: 



(30) 



' ^-4(l + -|)fM<-^|- (31) 



At that point, one has: 



2K 3 V if / V 2 



and the magnetization jmnp at the transition is: 

From this expression we immediately see that the condition of applicability of the Landau-Ginzburg expansion is 
J ^ {K -\- p), that is equivalent to consider a very stiff lattice or a very high pressure. Using the equation (|14|l . it is 
also possible to obtain an expression of the lattice parameter as: 

p J / 1\ p 3J2 

^=-f +xr-2j— f- 2.i.(ir+p) (^4) 

At the transition the jump of the lattice parameter is thus: 

= —6M = , ^. (35) 

if 27r2if(if+p) ^ ' 

The magnetic field at the transition is given by Eq. I|15() as: 

he p 3J2 



J is: %TT'^K{K+py 



(36) 



and obviously we have hc^ < he < hc^ . As we have already mentioned, in the intermediate regime, hysteresis can 
be observed as there are two minimas, one of them being metastable. This is shown on Fig. ^ where the energy for 
dE/dM = becomes superior to the energy of the fully saturated state ai h> he- At that point the thermodynamic 
minimum becomes Ai = 1/2, and the branch with dE/dM = becomes metastable. When the field is increased 
sufficiently rapidly, the magnetization may remain on the metastable branch until the field reaches /ic> • When the 
field is decreased from hc^ , at the solution with Ai =1/2 becomes metastable, and for a sufficiently rapid decrease 
of the magnetic field, one can have a magnetization Ai = 1/2 until the magnetic field reaches /ic<- This is the origin 
of the hysteresis cycle in the magnetization. In Fig. 0] is shown a plot of the hysteresis cycle of the magnetization for 
T = 0. A similar hysteresis cycle exists for the distortion and is observed in Fig. 3 of Ref. |3^ albeit in the vicinity of 
H = ■ In agreement with our picture, when the field is decreased, the system remains on the metastable distorted 
state, and when the distorted state becomes unstable the lattice parameter jumps to its value in the undistorted state. 
When the field is increased, the system remains in the metastable undistorted state until it becomes unstable. 

Returning to our model, we note that the application of pressure is pushing the first order transition to higher 
value of the applied magnetic field. This can be understood as a consequence of the fact that in the distorted state 
S < and \6\ is higher than in the undistorted state. As a result, the stability of the distorted state is increased by 
the application of pressure. Also, the hysteresis cycle becomes narrower as pressure is increased. An analogous effect 
is observed in ferromagnets^i, where application of pressure reestablishes the second order Curie transition. It might 
be interesting to investigate the effect of pressure on the transition studied in Ref. |3^ to see if the behavior of ii^ as 
a function of pressure is in agreement with our prediction. 
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FIG. 4: Hysteresis cycle for T = and K — J = 1, P = 0. The thermodynamic transition happens for h — he — 1.03J, 
however the distorted state is metastable for he < h < hc^ — 1.045J and the undistorted state is metastable in the region 
ftc< = J < h < he. As the field h is increased from to ^c> , the system remains on the metastable branch, and the 
magnetization jumps to 1/2 at ft = /!.c> • As the magnetic field is decreased, the magnetization drops from 1/2 at ft = ftc<. 
This produces an hysteresis cycle in the magnetization. A similar hysteresis cycle is expected for the lattice distortion. 



B. Behavior of susceptibilities in the distorted phase 

Now that we have described the first order transition, we would like to discuss the effect of the spin-phonon coupling 
on the susceptibilities of the system in the distorted phase. The following susceptibilities are of interest to us: 
the compressibility: 



(37) 

h 



the magnetic susceptibility: 



fdM\ 



and the parameter: 

'd5\ _ fdM 



measuring the Joule effect and its reciprocal, the Villari effedi^Si^. The two equivalent definitions of A are direct 
consequences of the definitions M ~ — (^), 6 = The equation of state results from the Eqs. H14I15|1 . 

Differentiating the second equation w.r.t. p, we find: 

fdM-^ ^ tan(^A/) f dS"^ ^^^^ 



\ dp J 7r(l — (5) \dp 



This equation indicates that applying a pressure for fixed magnetic field induces a change of magnetization (Villari 
effect), related to the compressibility. Differentiating the first equation w.r.t. p, instead gives: 
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Using Eq. (|40|l . we finally obtain the following expression for the compressibility: 



1 

K 



K — tan(7rM) sin(7rM) 



(42) 



We see that the compressibility of the system (which is in one dimension the inverse of the elastic constant) is enhanced 
by the interaction with the spin excitations. The physical reason is that when the lattice spacing is diminished, the 
magnetic energy is increased, which compensates for the loss of elastic energy. Expressed in terms of (5, using Eq. H15|l . 
the elastic constant reads: 



71- (1 - 5Y^J-^{l~5Y 



(43) 



The behavior of k is represented on Fig. [3 It presents a minimum for h = he- The jump of the elastic constant from 
a low value to its original value as the magnetic field is increased, implies a minimum in elastic constant as a function 
of magnetic field. At zero temperature such minimum appears more like a dip. Experimentally, the presence of such 
a dip near the saturation field is known as "A£"' effect in ferromagnetic materials^. Using Eqs. H36|l and H34|) . we 
obtain that the minimum value of the elastic constant reads for J <ti K: 



Kr — K 



1 



29 



J 



3 16 \tt{K +p) 



(44) 



showing that the elastic constant can be reduced to 1/3 of its value in the fully polarized phase. Applying pressure 
tends to diminish further the value of Kc- The behavior of k as a function of an applied magnetic field has been 
recently measured in a spin 1/2 Heisenberg chain material, the coordination polymer CuCCP'*''', via sound velocity 
measurements. These experiments reveal a pronounced minimum at the saturation field in qualitative agreement with 
the behavior represented on fig. |5| We will not attempt a quantitative comparison because the proper model to use 
to study the anomaly in the elastic constant of CuCCP is the Heisenberg model coupled to acoustic phonons and the 
thermodynamics of this model, even in the mean field approximation requires a much more sophisticated approach 
than the one of the present paperi^SiS^. 
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FIG. 5: The effective elastic constant as a function of the magnetic field. Below the saturation field, he, the elastic constant is 
very sharply reduced. For h = ha, the effective elastic constant jumps to K. The curve is drawn for J = K — 1. 
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The most characteristic signature of magnetostriction is the change of magnetization under apphed pressure. It is 
measured by: 



j^_(dM\ _ tan(7rM) _ / dS 

^ \~9p)h ^ T^Kil -5)~ Jtan(7rAf) sin(7rAf) ~ ^ \dh 



(45) 



The behavior of A is reported on FigEI The apphcation of pressure leads to a diminution of the magnetization, and 
the apphcation of a magnetic field in an increase of the lattice parameter, which indicates that the system has a 
positive magnetostriction^^. The variation of the magnetization A is an interesting quantity to measure, as well as 
(^), because they are non-zero only in the phase that is not fully polarized, and their equality is the very signature 
of magnetostriction. 

By considering the derivative of the equations of state H14I15|I w.r.t. h, we derive the magnetic susceptibility. 
Differentiating the first equation gives: 

d5\ J , , fdM\ , , 



Differentiating the second equation l(T^ then yields for the susceptibility: 



X - {dMIdh)^ ^ — , , . (47) 

ttJ [(1 — o) cos(7rM) — ^ sm (ttM)] 

The magnetic susceptibility can be rewritten as: 

X = Xo(/», J(1-'5))-, (48) 

K 

where Xo{h, J) is the magnetic susceptibility of an XX spin 1/2 chain with exchange constant J in a magnetic field 
h. We thus see that the susceptibility is strongly enhanced with respect to the case of a perfectly rigid lattice. A plot 
of X is reported in Fig|Bl The susceptibility has a maximum for h = he, and then falls to zero. Using Eqs. (|34|1 and 
we find that the maximum susceptibility is: 

_ K i 

Xc — ~i (.49] 

J I 189 / J 



32 i^iriK+p) 

and thus the maximum value of the susceptibility is proportional to the clastic modulus in the fully polarized phase. 
We also note that applying pressure results in an increase of the maximal susceptibility. With the expressions of A 
(Eq. 14011 and x l|47|l . we see that we have the following relation: 

{I- 5)- = -. 50 
X K 



IV. POSITIVE TEMPERATURE 



At fixed magnetic field and fixed pressure, we can vary the temperature and the numerical analysis of the free 
energy as a function of temperature shows that above a certain temperature, the metastable minimum disappears, 
and one is left with a single minimum. This behavior is illustrated for P — and = J on Fig. [S] As a result, 
above a certain temperature, the first order transition as a function of the magnetic field disappears and is replaced 
by a crossover. This is the analog of the critical point in the Van dor Waals mean field theory of the liquid-gas 
transition. In the present problem, the role of the pressure is played by the magnetic field. Applying the magnetic 
field is analogous to applying a pressure to the Van der Waals fluid to liquefy it. It is important to point out that the 
temperature at which the second minimum as a function of 5 disappears is a function of the applied magnetic field. 
The second minimum disappears at lower temperature when its energy is noticeably higher that the energy of the 
stable minimum. As a result, the width of the hysteresis cycle decreases as a function of temperature. A second effect 
of temperature is to blur the distinction between the fully magnetized phase and the partially magnetized phase. This 
effect can be seen on Fig. |H1 As temperature is increased, the stable minimum which corresponds to (5 = is moved 
to (5 < 0. Thus, as the magnetic field is decreased for T > 0, a reduction of the magnetization is observed even before 
the magnetization jumps, in contrast with the behavior for T = 0. 
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FIG. 6: The magnetic susceptibility as a function of the magnetic field. It diverges at the critical field where the transition 
takes place. The plot is done for J = K = 1. 
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FIG. 8: The evolution of the free energy as a function of temperature for K — J = 1 and h = 1.04. At zero temperature, the 
double minima arc neatly visible. The absolute minirrmm corresponds to 5 = indicating a non-magnetized phase, without 
deformation. For not too large positive temperature, the two minimum are still visible, but the barrier between them becomes 
smaller. Also, one notices that the minimum of the free energy corresponds to a solution with 5 < 0, indicating a contraction 
of the lattice as T increases, i.e. a negative dilatation coefficient. As temperature is further increased, i.e. for T = 0.03J, there 
is only a single minimum indicating the disappearance of the first order transition. 



A. Ginzburg-Landau expansion at finite T 

In this section, we turn to the derivation of the Ginzburg-Landau expansion at finite temperature. Such expansion 
is derived from the free energy: 



F{T, 5, h) = - i I g ln(l + 
in the vicinity of 5 = 0. Performing the Taylor expansion in the logarithmic function, we find: 



(51) 



F{T,S,h) 



dk 



e{k) 



^ j 2w e/^WA)-'') -h 1 



K- (3 



The solution with 6 = becomes unstable when: 



K{T) = K-f3 [ 



dk 



dk 



eikf 



cosh 



sinh 


fWfc)- 


h) 


5^ 


cosh^ 


■f(e(fc)- 


-h) 


24 



+ 0- 



cosh[j3{e{k) - h)] - 


2 (5* 


cosh^ 




192 



6(fc)2 



Stt 



cosh 



Ue{k)-h) 



< 0. 



At low temperature, we can estimate K{T) as: 



K{T) = K- - 



1 



/i2(j2 + 2/;/ 



(J2 - /l2)2 



(52) 



(53) 



(54) 
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Thus, the thermal correction can be neglected if: 



6 



2 

ft2 



(J2 - /l2)2 



« 1, 



(55) 



This leads for J ~ ft, to: 



T^^\J-h\ 



(56) 



Except for a small region of size of order tt^/ST around J, the thermal correction to the stability of the undistorted 
solution is negligible for low temperatures. This leads us to expect that the first order transition persists for non-zero 
but sufficiently low temperatures. We now show that there is a temperature r> above which the first-order transition 
as a function of the magnetic field disappears. Using a straightforward minoration cosh(2;) > 1 term in Eq. H53() . we 
find that: 



K^ — < K(T), 
8T ^ 



(57) 



and when: 



K > 



8T' 



(58) 



i.e. T > J2/(8JC), we have K{T) > and thus above a temperature T> = J'^/{8K), there is certainly no first order 
transition. As an aside, we note that the bound H58|l can be improved when h > J. Then, cosh[(e(fc) — h)/ (2T)] < 
cosh[(J- /i)/(2r)] and 



K 



STcosh^ {^] 



< K{T), 



(59) 



which shows that for low temperature, the solution with 6 = remains stable for high fields. We now proceed to show 
that for h < J an upper bound for K(T) can also be found in the regime of \h\ < J. Using this bound, we can show 
that the solution with (5 = is unstable at low but finite temperature on a soft lattice. This shows that the first order 
transition can be observed for T > in that case. To obtain the bound, we use kp defined by e{kp) = h and write: 



e(fc) 



- J(cos k — cos kp) = 2 J sin 



k 



kp\ . f k + kp 
sm 



(60) 



We have for fc > the obvious majoration: 



\e{k)-h\ <J\k-kp\ 



(61) 



Therefore: 



, 2 ^(k) — h\ o J\k — kp\ 

cosh^ ^ — - < cosh^ ' ' 



2T 



2T 



(62) 



and: 



K(T) < K - 



dk .P cos2 k 
^4cosh2^^^ 



(63) 



By a change of variables, this is rewritten as: 



K{T) < K 



2J riri^-kp) icos'^(kp + 



2T ^ 



dx- 



2T 



cosh^ : 



(64) 



In the limit Jkp ^ T and J(7r — kp) ^ T, it is possible to calculate approximately the upper bound as: 
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which shows that when K < J/(27r), the instabihty exists for low temperature T <C JMm{kp,Tr— kp)- The first order 
phase transition can therefore be observed at positive temperature at least when the lattice is sufficiently soft. 

We now turn to an estimation of the temperature at which the first order transition disappears valid in the case of 
a stiff lattice. First, we express K{T) using the density of states. We have: 

■■' de_ 

7rVJ2~~?4cosh2 (^) 

The divergence of the density of state for e ^ J is responsible for the instability. Therefore, it is legitimate, when 
\J — h\ <^ J to expand the integrand in Eq. for e ^ J. We find: 



K{T)^K- — ^ TT—TT (66) 



K{T)=K-l _ (67) 



TTy/2J{J -e) 4cosh^ 

Performing the change of variable e ^ J — { J — h)x, we can rewrite the integral in Eq. I|67fl as: 



73/2 fj-h^ 

K{T) = K-^==T C—A], (68) 



7r^2iJ-h) V T 
where: 

J^(^y)^r^ y (69) 

Jo \/^4cosh2 [f(l-x)] ^ ^ 

It is easy to see that J^{y 00) = 1 and T{y ^ 0) ~ y^^'^- As a result, at low temperature, T <C {J — h), we recover 
the zero temperature result, and at high temperature J 3> T ^ ( J — /i), we find that: 

73/2 

K(T)^K -7^, (70) 



and a T > Tc where: 



(71) 



the solution (5 = becomes stable, indicating that the first order transition disappears. We note that given the 
exponential temperature dependence of the spin-Peierls transition temperature^, the first order transition appears 
at a much higher temperature than the spin-Peierls transition in the case of a stiff lattice. 



B. Numerics 



By the minimization of the Gibbs free energy H10|l w.r.t. S, we have obtained the T dependence of the striction and 
the magnetization. 

We find that the discontinuity in S or M becomes smaller as T is increased. Also, we note there is a finite striction 
and a magnetization slightly below 1/2 in the high h phase. The reason for this is that at finite temperature, 
pseudofermion excitations are created which results in a lowering of the absolute value of the kinetic energy-like term 
in Eq. (|lll) . As a result, the absolute value of S is diminished in the partially magnetized phase and increased in the 
polarized phase. Above a certain temperature, the discontinuity is replaced by a large slope, indicating a continuous 
crossover from the fully polarized to the unpolarized state. This behavior is reminiscent of the liquid-gas transition 
in the Van der Waals fluid^^. For J = K = 1, we find the transition temperature to be of order 0.03 J in reasonable 
agreement with Eq. H71|l . 

C. Susceptibilities at positive temperature 

For positive temperatures, the equation of state is obtained from the equations: 

f dk e(k) 



1 -|_ e/3[(l-<5)e(fe)-''] 2 
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FIG. 9: The striction as a function of the magnetic field for various temperatures and K = J = 1. For T — O.Of, the 
discontinuity becomes smaller than for T — but it is still present. For T = 0.03 the discontinuity disappears, and is replaced 
by a rapid but continuous change of 6. For T = 0.1, the change becomes even more gradual. One can notice that 5 reaches 
zero at higher field when the temperature is increased. 
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FIG. 10: The magnetization as a function of the magnetic field for various temperatures and K = J — 1. For T = 0.01, the 
discontinuity is smaller, but it is still present. For T = O.OSdiscontinuity disappears, and is replaced by a rapid but continuous 
change of M. As temperature is further increased, the change of M becomes more gradual. One can notice that as T is 
increased, the saturation field becomes higher. 
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Taking the derivative of the first equation w.r.t. h yields: 



dh 



a f dk e(fc) 

_ 27r 4cosh2 |[(l--5)c(fc)-h] 



-A 



(74) 



The behavior of the Joule coefheient as a function of the magnetic field at finite temperature is reported in Figllll 
At very low temperatures it presents a cusp like behavior with a jump at the critical field he. As the temperature 
increases the cusp is replaced by a maximum for h = Hm (T) < he that becomes broader at increasing T. 
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FIG. 11: Derivative of the lattice parameter w.r.t. h and varying the temperature. The plot is done for J = K = 1. Below 
the transition temperature, the discontinuity is visible. Above the transition temperature, the discontinuity is replaced by a 
maximum. As the temperature increases, this maximum is lowered, and the region of anomaly is broadened. 



Taking the derivative of the second equation w.r.t. h then gives: 



( fl f ^ e(fc) 

/9M\ _ f dk 1 V J 27r4cosh2|[(l-5)e(fc)-;i] 



2 



2n 4cosh2 | [(l-5)c(fc)-h] 



\9hjj.p J 2^Ticosh^[{l~ S)e{k)-h] R ~ (3 J ^ ^ 



The plot of the magnetic susceptibility at finite temperature is reported in Fig. 1121 The singularity below Tc is changed 
into a maximum above T^. 

Similarly, taking the derivative of the first equation w.r.t. p yields: 

-1 _ /d5\ _ -1 

~'[d-p)^,,- K~f3j^--<^)^——' ^^^^ 

•' 2-7r 4cosh'^ f [(1— (5)e(fe)— /i] 

At finite temperature the effective elastic constant is thus given by: 



dk e{kf 
2^4cosh2f[(l-(5)e(fc) - h] 



«^^A'_^ / (77) 



As shown in Fig. E| the dip in the elastic constant at T = (cf. Fig.[SJ) as a function if the magnetic field is replaced 
by a minimum as the temperature T > Tc- The the T dependence of the elastic constant is very reminiscent of 
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FIG. 12: The magnetic susceptibility vs. h for varying T. At low T = 0.01, the magnetic susceptibility possesses a discontinuity. 
It is zero in the fully polarized phase, and becomes very large in the partially polarized state as the Fermi wavevector of magnetic 
excitations is nearing the band edge when magnetization is approaching saturation. For T = 0.03, the discontinuity is replaced 
by a peak. As temperature is further increased, the height of the peak diminishes, while its breadth increases. 



the results reported in Ref. |4J| with the dip becoming broader and less pronounced as the temperature is increased. 
However, in this reference, the behavior of the elastic constant was only discussed in a phenomenological way, by 
introducing an interaction proportional to the striction and the square of the magnetization, i.e. a term —S{{S^))'^. 
A more quantitative approach could be developed by treating the Hamiltonian by a TEA mean-field theory. 
Finally, taking the derivative of the second equation with respect to p allows one to check that Eq. (|74|) is recovered. 

D. Other thermodynamic quantities: Dilatation coefficient and specific heat 

Besides the susceptibilities considered in the zero temperature case, and recalculated at positive T, there are some 
other quantities that can be defined for finite temperature, namely the dilatation coefficient: 



the specific heat: 



where Ti. is the enthalpy, and the quantity: 



The internal energy is given by: 



d5 
df 



m 

df 



p.h 



(78) 



(79) 



1. Specific heat 



TJIT A M ^ f FiT,S,h) 



dk [{1 - S)e{k) - h] 
2^1 + e'3[(i-'5)^W-'»] 



(80) 
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FIG. 13: The effective elastic constant as a function of the magnetic field and different temperatures. The curve is drawn for 
J = K = 1. We notice that in the fully polarized state at low temperature, the elastic constant quickly returns to its value in 
the absence of interaction with magnetic excitations. A jump of elastic constant is observed as the magnetic field is crossing 
the saturation value. At higher temperature, the jump is replaced by a minimum in the elastic constant. The minimal value 
of the elastic constant is an increasing function of T. We also note that the width of the region where the elastic constant 
decreases is larger as T increases. 



The specific heat is given by the derivative of the internal energy, Cv = = — ^ (^)' El^pli*^itly, the specific 
heat at constant volume is given by: 



1 r dk [e{k) - hf 
J 2^ Acosh" §[e{k) - h]. 



To calculate the specific heat at constant pressure, we derive the enthalpy 7i (i, 5, h,p) = U (T, S, h) +pS and evaluate 
Cp{T,p) as: 



r dk e{k)[{1^5)eik)-h] 

^ I fdk [{l-5)e{k)-hr 1 1,^ 2-4cosh^( '->jr-"7 ^ 

-T^ J 2. 4 cosh- f [(1 - mk) -h]^T^K-l.J t ) ' 

where the second part comes from the variation of S with temperature. A plot of the specific heat is given on Fig. 1141 
for different temperatures. At very low-T the specific heat presents a singularity that is replaced by a maximum at 
finite temperature. 



2. Dilatation coefficient 



We use the definition of the dilatation coefficient H78|l . and Eq. lll|l to obtain: 

(83) 



r dk e(k)[{l-5)e(k)-h] 
J 2^4cosh^( <^-^>J<'°'-'' ) 



J 27T 



__e{kf__ 

4cosh^( '^-^'j],'°'-'- ) 
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FIG. 14: The specific heat for T = 0.01,0.03,0.05 and P = 0. We have K = J = 1. The specific heat becomes very large 
near the transition. Above the transition temperature, a peak can be observed in the specific heat. The height of this peak 
diminishes as T increases, whereas its breadth increases. 



We note from Eq. H83|) that the usual relation between Cp, Cy, the compressibility k and a (Eq. (16,9) of Ref . Is^l is 
satisfied. The variation of the dilatation coefficient with magnetic field is shown in Fig^l As we see, when the first 
order transition is present, such coefficient vanishes in the fully polarized phase and becomes negative in the partially 
polarized phase. Above the temperature Tc, the dilatation coefficient varies continuously, but its behavior is non- 
monotonous when the magnetic field is close to the value that induced the first-order transition at low temperature. 



V. CONCLUSIONS 



In this paper, we have discussed magnetostriction in an array of XX spin chains coupled to three dimensional 
acoustic phonons. We have shown that for low temperature, a first order transition would be obtained as a function 
of the magnetic field, with a jump of both the magnetization and the lattice constant. For higher temperatures, the 
discontinuity is replaced by a crossover, but anomalies are still visible in the elastic constant, the specific heat, the 
dilatation coefficient and the magnetic susceptibility. As we have explained in the introduction, our results are also 
relevant to ladder systems under magnetic field since the latter can be mapped in the limit of strong rung coupling onto 
XXZ chain model under field. We stress again the important differences between magnetostriction and the generalized 
spin Peierls transition of Ref. l37l First, in the generalized spin Peierls transition, a static lattice deformation of 
wavevector q = 2-kM where M is the magnetization is expected. This would lead to extra reflections in an X-Ray 
scattering experiment. On the contrary, in magnetostriction, no superlattice is formed, and no new reflections are 
expected. Second, in the generalized spin Peierls transition, a spin gap is formed, leading to an exponential suppression 
of the specific heat and magnetic susceptibility at low temperature. In the case of magnetostriction, no spin gap is 
formed, leading to a T linear specific heat and a finite susceptibility at low temperature in the partially magnetized 
phase. Also, in the case of the magnetostriction, anomalies are expected in the elastic constant and the dilatation 
coefficient as a function of the magnetic field. Of course, in a real system, both optical and acoustic phonons are 
present, and one should expect a coexistence of magnetostriction and generalized spin-Peierls effects^. However, the 
first order transition under applied field is observable at higher temperature than the spin-Peierls transition so that the 
two effects can be separated in principle. It would be interesting to analyze the more general case of magnetostriction 
effects in a XXZ chain by using the TBA^iiSiSi^ or by using Quantum Monte Carlo method. We expect that the 
general conclusions concerning the existence of the magnetostriction effect and the possibility to separate it from 
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FIG. 15: Dilatation coefficient for T = 0.01,0.03,0.05 for P = and A" = J = 1 as a function of the magnetic field. The 
dilatation coefficient vanishes in the fully polarized phase, and is negative in the partially polarized phase. We notice that the 
dilatation coefficient decreases strongly as temperature is increased. The behavior of the dilatation coefficient is non-monotonous 
as a function of the magnetic field. 



the spin-Peierls transition will not be affected, and a more quantitative agreement with experiments will obtain. A 
more interesting extension of the present work would be to consider a two dimensional system of hard bosons coupled 
to phonons as this is relevant to the experiments on the TlCuCla material^ and the CuHpCl material^^. In two 
dimensions, hard bosons on a lattice can be mapped onto fermions in a gauge fielc&^. For a low fermion density, 
the corresponding gauge field is small. Therefore, at low magnetization, it can be expected that one can describe 
magnetostriction effect in a manner similar to the ID case. We will leave this for future work. 
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CALCULATION OF THE FIELD h^^ 

The Eq. ipn|l can be rewritten as: 

P{S) = [n^KS + pf - J2] (1 - 5f = 
The extrema of the function P are given by dP/{dS) ~ i.e.: 



giving a maximum for S — 5 + +, with: 

1 3p 
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K \ttK 
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= 0, 



(84) 
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and a minimum for S = 5^, with: 



(5_ = 



1 - 



3p 
K 
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ttK J 



(87) 



The critical field h^^ is obtained by solving simultaneously the equation: 



P{6-) = 

This is done by introducing the polynomial: 

P{5) = 52 + {-ip/K - 1)5 + 



hi 



K 



K 



J 



(89) 



and writing P{S) - hl^ = Qi6)PiS) + R{6). Since F((5_) = h^^ and P((5_) = 0, the equation lS^_reduced to 
R{S-) — 0. Since i? is a polynomial of degree 1, such an equation is trivial to solve. This leads to Eq. 
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